Measurements of critical current diffraction patterns 
in annular Josephson junctions 
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We report systematic measurements of the critical current versus magnetic field patterns of annular 
Josephson junctions in a wide magnetic field range. A modulation of the envelope of the pattern, 
which depends on the junction width, is observed. The data are compared with theory and good 
agreement is found. 
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Large area Josephson junctions are intriguing objects 
for performing experiments on nonlinear electrodynam- 
ics. In particular, the propagation of solitons, also called 
Josephson vortices or fluxons, has attracted a lot of at- 
tention and has been studied in detailB Recently large 
area Josephson junctions have been proposecLto be used 
as efficient radiation and particle detectors.Enl In such 
junctions, the spatial dependence of the phase difference 
between the superconducting electrodes is an important 
characteristic that determines the junction properties. 

In an annular Josephson junction, magnetic flux 
quanta threading one superconducting loop but not the 
other, are trapped anjd stored in the junction due to the 
fluxoid quantization.^ This property of the system offers 
the unique possibility to study fluxon_dynamics in the 
absence of collisions-. with boundaries In the annular 
junctions proposecoQ for radiation and particle detec- 
tion, trapped vortices are useful to suppress the critical 
current of the junction in order to allow for a stable bias 
point in the subgap region of the junction current- voltage 
characteristic. 

In this report, we present systematic measurements 
of the critical current of annular Josephson junctions 
in dependence on the externally applied in-plane mag- 
netic field. The critical current I c of a junction with- 
out trapped fluxons is at maximum when no magnetic 
fields are present. In the presence of a magnetic field 
this maximum superconducting current is reduced. Mag- 
netic fields can be due to the bias current applied to the 
junction (self-fields), due to flux trapped in the junction 
itself or its leads (Josephson or Abrikosov vortices, re- 
spectively), or they can be applied externally. The mod- 
ulation of the critical current with the external field is 
often called a critical current diffraction pattern. We in- 
vestigate these patterns for annular junctions of various 
dimensions in a wide range of magnetic fields. 

We present experimental data on five annular Joseph- 
son junctions with the same external radius r e — 50 /im 
but different inner radii ranging from 30 to 47 /im, 
see second and third column of Tab. |. Hence the width 
w = r e — Ti of the junctions is ranging from 3 to 20 ^m. 
The junction geometry is shown in Fig. All junc- 
tions have been prepared on the same chip using Hypres 



TABLE I. Geometrical parameters and fitted values of the 
measured annular Josephson junctions. 



junction # 


r-i [fim] 


5 — n /r e 


AH' /AH 2r e /w 


5 


Ar [/im 


A 


47 


0.94 






0.96 


0.5 


B 


45 


0.9 


22.9 


20.0 


0.92 


0.5 


C 


42 


0.84 


13.2 


12.5 


0.88 


1.0 


D 


35 


0.7 


6.5 


6.7 


0.72 


0.6 


E 


30 


0.6 


4.8 


5.0 


0.62 


0.5 


junction # Ho [Oe] 


Ho /Ho 


A [nm] 








A 


0.319 


0.65 


208 








B 


0.346 


0.703 


193 








C 


0.321 


0.646 


210 








D 


0.405 


0.821 


165 








E 


0.376 


0.765 


177 









technology! with a nominal critical current density of 
j c = 100 A/cm 2 . Accordingly, the Josephson length is 
approximately 30 /im at 4.2 K. 

In Fig. H the critical current diffraction patterns of 
the two junctions B and D, being representative for 
the set of measured samples, are shown. Obviously, a 
strong dependence of the pattern on the junction width 
is observed. As expected, the critical current of the 
junction at zero field scales with the junction size as 
I c = j c 27r(r 2 — rf). Measuring the diffraction patterns in 
a wide range of magnetic field, two characteristic modu- 
lation scales of the critical current are observed. The pat- 
tern having a small magnetic field period AH has an en- 
velope of the larger period AH' which depends strongly 
on the width of the junction (compare Figs. ||a and b). 

The observed critical current diffraction patterns can 
be qualitatively understood in the following way. The 
modulation of the period AH is due to the penetration of 
magnetic flux in the direction perpendicular to the exter- 
nal magnetic field. This period is inversely proportional 
to the diameter of the junction: AH oc l/(2r e ). This 
is analogous to the standard case, where AH is propor- 
tional to the reciprocal junctions length in the direction 
perpendicular to the magnetic field.cl The minima of the 
modulation of the period AH' occur, when the magnetic 
flux penetrates the junction strongly also along the width 
of the junction. Therefore, the period AH' of the second 
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FIG. 1. Schematic of the Annular Josephson junction bi- 
ased in the Lyngby-geometryfl (a) and its dimensions (b). 



modulation is proportional to 1/w. By calculating the 
ratio 



AH' 
~AH 



2re 
w 



(1) 



for the different junctions, this simple prediction can be 
quantitatively compared with experiment. As can be 
seen from the forth and fifth column of Tab. |, Eq. (Q) is 
quite accurately fulfilled for our junctions!] 

The described effect can be illustrated by plotting the 
supercurrent density j s at different magnetic magnetic 
fields versus the junction coordinates. At the magnetic 
field H — 3.25 Oe < AH' , approximately two and a half 
flux quanta penetrate into the junction cross section 2r e , 
as shown in the inset I of Fig. At the larger field 
H = 12.2 Oe > AH' , more than one flux quantum pene- 
trates the width cross section of the junction (see inset II 
of Fig. ||b). Thus, after each period AH', one additional 
flux quantum has penetrated the width of the junction. 
We note here, that the spatial distribution of the supet 
current density could also be measured in experiment JH3 

Several approaches to calculate the critical current 
diffraction patterns of aaruijar Josephson junctions have 
been published earlier. tT£3 Mainly, two different cases 
have been considered, the long annular Josephson junc- 
tion with a circumference 27rr larger than the Joseph- 
son length Xj and-the small annular Josephson junction 
where 2nr < Xjf3 where r = (rj + r e )/2 is the mean 
radius of the junction. 

The most complete theoretical description of the crit- 
ical current diffraction pattern I C {H) of small annular 
junctions of arbitrary width and number of trapped flux- 
ons n is presented by Nappi in Rcf. 13. The dependence 
of the critical current I c on the magnetic field H is given 
by the formula 



Ir = h 



T~P l xJn ( x lk 1 dx 




H (Oe) 

FIG. 2. Critical current diffraction patterns of (a) junction 
B and (b) junction D at 4.2 K. Dots are experimental data, 
the solid line is theory according to Eq. (^). For better vis- 
ibility the low current region is also plotted on an enlarged 
scale. The two field modulation periods AH and AH' are 
indicated in each plot. The insets / and II of plot b) display 
the supercurrent distribution in junction D at the magnetic 
fields indicated by arrows; light (dark) regions correspond to 
current in positive (negative) direction. 



(2) 
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where J„ is n-th Bessel function of integer order, 6 = 
Ti I r e is the ratio of inner radius r, to outer radius r e and 
Iq is the maximum superconducting current at zero field. 
The field 

Ho = $o/(27rr eW) A) (3) 

is the characteristic magnetic field; $o is the flux quan- 
tum, /i is the vacuum permeability and A is the effective 
magnetic thickness. lj 

For n = 0, the two extreme cases 5 — ► 1 (see Ref. [fl]) 
and 6 — ► (Ref. ||) of Eq. (||) have been discussed in 
the literature. The predictions of Eq. (||) have also been 
compared topfixpjeriments in a relatively small magnetic 
field range D'EJ'Ej To our knowledge, there has been no 
systematic comparison of the theory with experimental 
data for different junction width in a large field range. 
Our intention here is to perform such a comparison. 

In Figure ^, our experimental data are fitted to Eq. (g). 
In the fitting procedure the values of both H and 5 are 
determined. Subsequently, the quantities acquired from 
the fits are labeled by a tilde (Ho, 8). For the fit, the 
initial value of 5 is calculated from the designed geome- 
try of the junction; the initial Ho is calculated according 
to Eq. (|}|) assuming the reasonable value of 200 nm for 
the magnetic thickness A. Then, the best fit is found by 
iteratively adjusting Hq and <5. The value of H$ predomi- 
nantly determines the small period of the critical current 
modulation AH, whereas S determines the large mod- 
ulation scale AH' . This fact is in agreement with the 
qualitative discussion above. As can be seen from Fig. ||, 
excellent agreement between theory and experiment is 
found. The parameters 5 and Hq determined from the 
best fits to the data of junctions A to E are quoted in 
Tab. | 

Comparing the values of S and S in Tab. we find that 
S < 5 for all junctions. This small but systematic devia- 
tion can be explained by assuming a symmetric deviation 
Ar of the junction radii from their designed dimensions 
(e.g., due to the photolithographic procedure during the 
preparation). Using this assumption, 8 can be expressed 
as 

6=^. (4) 

From the fits we find that Ar varies between 0.5 and 
1.0 /im (see Tab. |). This size correction can also be ex- 
plained as due to a slight over-etching of the trilayers 
during sample fabrication, that results in a small reduc- 
tion of the sample size. The obtained Ax values agree 
with the size tolerance quoted by HypresE 

According to theory, the quantity Hq does only depend 
on the outer junction radius r e and hence should be iden- 
tical for all junctions measured. Instead, we find values of 
Hq from the fits that slightly vary from junction to junc- 
tion, see Tab. [j. Using Eq. (||) the magnetic thickness 
A can be calculated from Hq. The values of A obtained 
for each junction are quoted in the last column of Tab. Q. 



The average magnetic thickness is A = 191±18 nm, which 
is in good agreement with the value of A ss 2A^ in the 
thick film limit, yielding a London penetration depth of 
Xl ~ 95 nm. 

The scatter observed in Hq (or, equivalently, in A) may 
be due to a small number m of flux quanta threading the 
holes of both junction electrodes simultaneously. The 
critical current diffraction patterns for different values m 
are very similar in their qualitative features, but may dif- 
fer quantitatively. Preliminary experimental results show 
that, upon cooling the junction from the normal to the 
superconducting state in a small residual magnetic field a 
large number of times and measuring the resulting critical 
current versus magnetic field, slightly different diffraction 
patterns depending on the value of m are observed. In 
such measurements we have only observed three differ- 
ent diffraction patterns, despite repeating the described 
procedure a large number of times. This strongly sug- 
gests that this effect is due to magnetic flux threading 
the junction loop perpendicular to the substrate. 

At small fields, we observe a systematic deviation of 
the calculated patterns from the experimental ones. In 
particular, the first minimum of the critical current ap- 
pears at larger field values than predicted by the the- 
ory. Moreover, the critical current at the first minimum 
does not fall to zero. Both facts are to be expected for 
junctions that are not really small in comparison with 
A,/. Indeed, the dimensions of our junctions are slightly 
larger than Aj. This leads to a inhomogeneous penetra- 
tion of the magnetic field into the junction at low fields, 
resulting in an increase of the field value H at which the 
the first minimum of the pattern is observed. The anal- 
ogous effect-is observed in conventional long Josephson 
j unctions E'ED 

At higher temperatures, the Josephson length Aj 
increasesB and, hence, the effective size of the junction de- 
creases. In the inset of Fig. [|, the calculated normalized 
external junction radius r e /Xj is plotted versus temper- 
ature, taking into account the temperature dependence 
of both the critical current. density j c {T) and the London 
penetration depth \l(T)H At T > 7.8 K the normalized 
radius drops below unity. Therefore, at higher tempera- 
tures a better agreement between experimental data and 
theory can be expected at low fields. This is illustrated 
in Fig. [§[ where the experimental critical current diffrac- 
tion pattern of junction B is plotted together with a fit 
for the temperatures T = 4.0, 7.0, 8.5 K. The fit is made 
keeping 5 constant for all T and adjusting Hq. At ele- 
vated temperatures, both the position of the first mini- 
mum and the modulation depth of the critical current at 
small fields show better agreement with the theoretical 
prediction. 

We have also measured junctions with a single fluxon 
trapped in the junction barrier (n — I). As an example, 
the critical current diffraction pattern of junction B for 
n = f at 4.2 K is shown in Fig. [|. Taking the same fitting 
parameters as for the case of no trapped fluxon (n = 
0), we find as good agreement between the theory and 
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FIG. 3. Critical current diffraction patterns of junction B 
at temperatures between 4K and 8.5K. Dots are experimen- 
tal data, solid line is theory. In the inset the calculated nor- 
malized external junction radius r e /Xj is plotted versus tem- 
perature. 
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FIG. 4. Diffraction pattern of junction B with one trapped 
fluxon at 4.2 K. Dots are experimental data, solid line is the- 
ory. The inset displays the supercurrent distribution in the 
junction at the magnetic field indicated by the arrow. 



the experimental data as before. The slight differences 
between the fit and the experimental data are, again, due 
to the dimensions (r e > Aj) of the junction. In the inset 
of Fig. ^ the supercurrent distribution in junction B at 
H = 1.71 Oe calculated according to Eq. (||) is shown. 
Obviously at this field a number of vortex anti-vortex 
pairs have penetrated into the junction but the width 
of the junction is not fully penetrated (compare Fig. ||b, 
inset II). The symmetry of the current distribution in 
the junction is broken due to the presence of the trapped 
vortex. Similar current distributions in the presence pfj 
trapped vortices have also been observed in experiment .E2I 
It is worth to point out that good agreement between 
theory and experiment in the large field range is found 
for junctions of a diameter substantially larger than Aj. 
At low fields the theoryEj describes well the experiments 
with r e < A j, as confirmed by our measurements at 
higher temperatures. Thus, the magnetic properties of 
the junction are determined rather by the junction radius 
than by the junction circumference, as already pointed 
out in Ref. 

In summary, we have systematically measured the crit- 
ical current diffraction patterns of a number of annular 
junctions of different width, with and without trapped 
fluxons, in a wide magnetic field range and at different 
temperatures. The experimental data show a pronounced 
width dependence that is explained accurately using the 
existing theory. In particular, a modulation of the enve- 
lope of the critical current diffraction pattern is observed 
for junctions of large width. The period of this modula- 
tion depends very sensitively on the normalized junction 
size described by the parameter 5. The method of our 
data analysis is accurate enough to detect a small re- 



duction of the size of the junction due to the fabrication 
process. 
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